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Abstract 



Use is made of Onsager's hydro dynamic equation to derive the vibration 
spectrum of the vortex lattice in d-wave superconductor. In particular the 
rhombic lattice (i.e. the 45° tilted square lattice) is found to be stable for 
B > H cr (t). Here H cr (t) denotes the critical field at which the vortex lattice 
transition takes place. 
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Typeset using REVT^X 



Since the discovery of the triangular vortex lattice in type II superconductors by 
Abrikosov Q and others fHH, the vibrational modes of the vortex lattice have been studied 
by a number of people H. In these works the stability of the triangular vortex lattice in an 
s-wave superconductor is established. 

The discovery of hole-doped high T c cuprate superconductivity by Bednorz and Miiller 
H in 1986 and the recent realization |J that <i-wave superconductivity is involved may give 
a new twist on the whole subject. 

We have shown earlier that the rhombic vortex lattice (or the 45° tilted square lattice) is 
stable in the vicinity of the upper critical field in a (i-wave superconductor in a magnetic field 
parallel to the c axis |7J . More recently we have shown || the vortex lattice transition from 
the triangular lattice to the square lattice takes place at a small magnetic field H cr (t) ~ 
K^H&^t), which implies that the vortex lattice should be rhombic in the overwhelming 
region in the B-T phase diagram (see Fig. 1). Here k is the Ginzburg-Landau parameter 
and we have H cr about a few Tesla in YBCO and Bi2212 at low temperatures. Indeed this 
critical field H cr (t) is consistent with the observation of rhombic vortex lattice by SANS 



and by STM imaging |10[ in monocrystals of YBCO at 3 Tesla. 

The object of this paper is to study the vibrational spectrum of the 45° tilted square 
vortex lattice in a magnetic field. First following Fetter et al. M, we study Onsager's and 



Landau's hydrodynamic equation ]TT],12] for a vortex lattice. Basically we assume that 



the vortex moves with the local velocity generated by other vortices. We find that the 
square vortex lattice is unstable when B < H cr (t) but becomes stable for B > H cr (t). 
The vibrational spectrum are determined in the whole Brillouin zone. Second we analyze 
the vibration spectrum within the time dependent Ginzburg-Landau equation |TJJ with the 



Aranov-Hikami-Larkin term ||15|| . In this limit we have now a set of damped oscillation 
modes rather than oscillation modes. However the stability condition of the square vortex 
lattice is the same as in the analysis using Landau's hydrodynamic equation as expected. For 
example this new vibration spectrum will be crucial in determining the melting transition 
line where the vortex lattice melts into a vortex liquid. 

I. Vibration modes in the square lattice (hydrodynamic limit) 

As is well known the Landau's and Onsager's hydrodynamic equation applies when motion 
of vortices involves no energy dissipation. Unfortunately this condition is never realized for 
vortices in a d-wave superconductor since there are the low energy extended states attached 



to every vortex and they certainly dissipate energy whenever the vortex is in motion ||T6 
Nevertheless it is of great interest to study this idealized limit. 

The extended Ginzburg-Landau free energy for the vortex state with ^ < (I < A in 
o?-wave superconductors reduces to M 

tt=^^£V(r,-ro), (1) 

K L 

0(r) = - e a 2 £V— + K [-) , (2) 
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where e = 31C(5)(- lnt)/196C(3) 2 ~ 0.114(-lnt), t = T c /T, a 2 = f (In 2 + §) = 5.454 •• {]; 
and £, d, and A are the coherence length, intervortex distance and the magnetic penetration 
depth respectively. Here 0(r) describes the interaction energy between two vortices; the first 
term is the core interaction (new to d-w&ve superconductors) while the second term is the 
usual magnetic interaction. Following Onsager |12| the equation for the L-th vortex is 



2ty 



[z x V) ^lu ■ V) (t>{v L , - V L ) 



(3) 



L'+L 



where R = h/2m = ir/m (in quantum unit), yl is the position of the L-the vortex, is 
the small displacement from the equilibrium position r° L (namely tl = ul + r°), V denotes 
the differentiation with respect to r^. We have set vllv = - ul>, ylv — y l — rjy- In 
deriving Eq.(^|) the Hamiltonian H(tl) in Onsager's theory is replaced by the free energy 
£1(tl) obtained in Ref. ||. In Ref. ||, the square vortex lattice tilted 45° was studied. In 
the following analysis, however, it is more convenient to use the coordinates system rotated 
by 7r/4 around the z axis, so that vortex lattice axes coincide with the x- and y-axes. This 
change makes the sign in front of the core interaction term opposite from the one in Ref. 
||. The equilibrium positions of the square vortex lattice become r° L = d(m,n) (m,n G Z), 
where d = yJfo/B (0 O is the flux quantum). 

Now introducing a plane wave representation as u L = sexp(i(q • r° L — out)), Eq.(^) is 
rewritten as 



— iuos x = as x + /3s y , 
-iusy = is x - as y , 



(4) 



where 



a 



2ttA : 



E'(i-^; 



r 



^ 9 . 9 sin 69 1 
, K 2 ( -f ) - 20ea 2 K 2 e 4 A 2 — ^ 
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r. 



rz, • 



T) 
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4ttA 2 



47rA : 



= — A2 -r + 40ea 2 /t 4 A z s — 

rf V A / r? 



E'(i-^)^ (y) 



(5) 
(6) 
(7) 



where i^o^) an d K 2 (x) are the modified Bessel functions. The vibration frequency u> is 
then given by 
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-a A - /3 7 = r]' z -a 2 - £ 2 



(8) 



Let us consider the long wave length limit (q = |q| <C 1) of the vibration spectrum. Then 
we obtain 



1 Note that a 2 is related to the constant a\ given in Ref. || as a 2 = 5ai/2ir. 



3 



— £ + ia 



V 



2 J2 



32ttA : 



16ttA 2 

where x is the angle q makes with the x axis, we have set 



2f4 



a(B) 

where \i = d/A(<C 1) and 
H CT (t) = c 



40ea 2 K / £ 



■fi 2 ~ 0.145025 



B 

H rr 



AnC- 1 



U V 10ea 2K 2 e(A8 + A 4 ) 
We have also used the notations 



0.48(-lnt) 



-1/2 H c2 (t) 



^i = d- 2 Y,'rlMr L /\) = - A +0{l) 

L ^ 

S 24 = d~^'rlK 2 (r L /X) = ^-^ + 0(1), 

L 



/I 4 /i 2 



L 'L P P 

A 8 = rf 4 ^ '_ C os(8fl L ) = 5.0306 
A 4 = d 4 yl cos(4fl L ) = 3. 1512 ■ ■ ■ , 



0(1), 



(9) 
(10) 

(11) 

(12) 

(13) 
(14) 
(15) 
(16) 
(17) 



where C = 0.10331. The quantities Si, X 24 , £^4 and C are quoted from || for the sake of 
the readers convenience. Using these, we have 



327T 

Rq 2 H 2 j 



32tt 



[2 + a(B)A 8 ) - (16nC - 2 - a(5)A 4 ) cos4x + 



(18) 



Finally, we obtain 

( 



uj 2 ~ 



(eB 
mc 



q 4 X 2 d 2 



( B\ CT ) 



0.01989 + 0.007257 
- f 0.03176 — 0.0045 IG 



B 



COS 4% 



(19) 



The condition to have real spectrum (uj 2 > 0) is B > H cr . We see immediately that the 
square vortex lattice is unstable for B < H cr for some \i while it becomes stable for all \ for 
B > H cr as it should be. Exactly at the transition point B = H cr uj 2 = at \ = 0, vr/2, tt 
and 37r/2. 

In the previous paper ||, we derived the critical field as H cr {t) ~ 0.52(— hit) -1 / 2 H c2 (t) / ' k 
by minimizing the free energy of the extended Landau Ginzburg theory. This and the result 
obtained here in Eq. ([12]) appear to be fully consistent. 
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II. Thermal fluctuation of the vortex lattice 

Let us consider thermal fluctuation of the square vortex lattice (B > H cr ) based on the 
Landau- Ginzburg free energy ([I]). This provides us with another way to study the quantities 
a, P and 7 introduced in the last section. 

The small deviation of the free energy AQ due to the lattice vibration is 



AQ 



im^ 2 



E ' E _ «o M ) < ^V(r ( L )( u Lu - u Qu ) 



(20) 



L H,v=x,y 



where denotes the //-component of the small deviation vector u L and (^(r) = d^d u (j){r). 
Introducing the Fourier transformation Ul^ = D q e * q r <v(q) (^ ne rea lity condition s M (— q) = 
s M (q)* is imposed), we have 



AQ 



4vr 2 n^ 2 

K 2 K 

4vr 2 ^ 2 

hZ 2 R 



—7 a 

a (3 



E(^+k+(q)| 2 + ^|s_(q) 



(21) 



where s±(q) denotes a suitable unitary transformation of s M (q) and 



cu± = i] ± J a 2 + £ 



K 



32nX 



16vr /16tt (2 + a(5)A 8 ) - (16ttC - 2 - a(5)A 4 ) cos 4;^ 



(-) 

\mcj 



1±1 (d" 



0.01989 + 0.007257 



B 



- 0.03176 + 0.004546 



B 



cos A\ 



(22) 



(23) 



The stability condition u;± > leads to 77 > y^ 2 + £ 2 or B > H cr . Note that this condition is 
exactly the same as the one in the previous section. Evaluating the path integral / Ds (\s x \ 2 + 
\s y \ 2 )e-P n / J Vs e _/3n , we obtain 



E(Mq)| 2 + l%(q)| 2 > = A:T 



2 — 

K K 



E 



1 



1 



(24) 



Studying this expression we will be able to study the vortex lattice melting from the square 



vortex lattice. Details on the vortex lattice melting will be reported in a separate paper [13 



III. Time dependent Ginzburg Landau equation approach 

In a (i-wave superconductor the vortex motion cannot be purely hydrodynamic, since the 
vortex motion always accompanied with the energy dissipation as seen from a number of 
experiments involving the flux flow resistance in high T c cuprate superconductors. Therefore 
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the vibration models for the vortex lattice should be described better by the time- dependent 
Ginzburg Landau equation [14], supplemented by the term which accounts for the Hall effect 
| T5| at least qualitatively. 

In the vicinity of T = T c the time dependent Ginzburg Landau (TDGL) equation is the 
same as that in s-wave superconductors |T^]. Neglecting the term which gives rise to the 
Hall effect we obtain 



-A(r) = iV 1 



8T r dt 



dA*(r) 



(25) 



where Q is the free energy defined in Ref. M. Assuming the usual product expression 
A(r) = ATJi /(r — r.j), we obtain 



71 a 2 1 

A 2 v- 

8TL 2 



df(r - r; 



dr.; 



ni/( r -^ 



^ dv,-; 



(26) 



which reduces to 



_ 16T c e(T) 

V,; — . „ iV n 



7i A 2 (T) Dri 



(27) 



The right hand side of Eq. (p7|) is evaluated as in the case of the hydrodynamic limit, and 
obtain two frequencies 



u ± = -iA Q ( V ± Je + a 2 ) 



where rj ± \/^ 2 + a 2 is calculated in Eq. ( |2^ ) and 

U((3)E F 



An 



7T 



3 Tr 



(28) 



(29) 



The stability of the vortex lattice requires rj > y/C, 2 + a' 2 , which is satisfied for B > H cr (t). 
Finally the gauge invariance of TDGL requires that Eq.(p5|) has to be replaced by 



where 



(l-„)_A(r) = V™ 



4T r d\nT r 



7T 9/i 

Then in the presence of the Hall effect Eq.(|2B]) should be replaced by 

a — i 

UJ± 



(30) 



(31) 



l + a 2 



A ( V ±Je + a 2 ). 



(32) 



Therefore in this general situation the vibration becomes damped oscillation. However, the 
stability condition is the same as before. 
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Concluding Remarks 

Limiting ourselves to the square vortex lattice, we have studied the oscillation mode of 
the vortex lattice first by the hydrodynamic equation and then by the time-dependent 
Ginzburg Landau equation. Both analysis indicates that the square vortex lattice is stable 
for B > H cr (t), while it becomes unstable for B < H cr (t). Of particular interest is that 
the fluctuation of the vortex lattice diverges linearly like (B — iJ cr (t)) -1 / 2 as B approaches 
the transition point. This should have a profound implication on the melting of the vortex 
lattice for example. The related problems will be discussed in future publications. 

Acknowledgment 

One of us (KM) would like to thank CREST enabled him to visit two weeks at ISSP, 
University of Tokyo. This work is supported by NSF under grant number DMR 95-31720. 



7 



REFERENCES 



[1] A. A. Abrikosov, Soviet Phys. JETP 5 1174 (1957). 
[2] J. Matricon, Phys. Lett 9 289 (1964). 

[3] W.H. Kleiner, L.M. Roth and S.H. Autler, Phys. Rev. 133 A1226 (1964); G. Eilen- 

berger, Z. Physik 180 32 (1964); D. Saint- James, G. Sarma and E.J. Thomas, "Type II 

superconductivity" (Pergamon, London 1968). 
[4] See for example, A.L. Fetter, P.C. Hohenberg and P. Pincus, Phys. Rev. 147 140 (1966). 
[5] J.G. Bednorz and K.A.Muller, Z. Phys. B 64 189 (1986). 
[6] See for example K. Maki and H. Won, J. Phys. I(France) 6 2317(1996). 
[7] H. Won and K. Maki, Europhys. Lett. 30 421 (1995); Phys. Rev. B53 5927 (1996). 
[8] J. Shiraishi, M. Kohmoto and K. Maki, Phys. Rev. B (in press); also in "Symmetry and 

Pairing in Superconductors" edited by M. Ausloos, (Kluwer, Dordrecht 1998). 
[9] B. Keimer, W.Y. Shih, RW. Lynn, F. Dogan and LA. Aksay, Phys. Rev. Lett. 73 3459 

(1994). 

[10] I. Maggio-Aprile, Ch. Renner, A. Erb, E. Walker and 0. Fischer, Phys. Rev. Lett. 75 

2754 (1995). 
[11] L. D. Landau, J. Phys. (USSR) 5 71 (1941). 
[12] L. Onsager, Nuvo Cimento 6 suppl II 249 (1949). 
[13] in preparation. 

[14] A. Schmid, Phys. Kond. Materie. 5 302 (19660; C. Caroli and K. Maki, Phys. Rev. 159 
306 (1967). 

[15] A.G. Aranov, S. Hikami and A.I. Larkin, Phys. Rev. B51 3880 (1995); H. Fukuyama, 
H. Ebisawa and T. Tsuzuki, Prog. Theor. Phys. 46 1028 (1971). 

[16] Y. Morita, M. Kohmoto and K. Maki, Phys. Rev. Lett. 78 4841 (1997); Phys. Rev. Lett. 
79 4514(c) (1997); Europhys. Lett. 70 207 (1997). 



8 



FIGURES 

Fig. 1 The B-T phsae diagram. 
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